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Abstract: The results from the first three years observation in Telescope Array are reviewed. The energy spectrum, mass 
composition, and anisotropy in arrival directions of ultra-high energy cosmic rays are discussed. The energy spectrum 
of cosmic rays with energies above 10 18 eV is determined using the data of the fluorescence detectors, surface detectors, 
and with a hybrid shower reconstruction technique using the both data. The ankle structure is clearly found at 10 18 7 eV. 
The second bending point is also identified, at which the flux of cosmic rays begins to decrease and the spectral index 
changes from E~ 2 ' 68 to E~ 4 ' 2 . The study of shower development using the fluorescence detector data shows that proton 
is dominant in cosmic rays above 10 18-2 eV. No apparent anisotropics in the arrival direction of cosmic rays are found. I 
also present further experimental extensions which are being developed at the Telescope Array site. 

Keywords: Ultra-high energy cosmic rays, Telescope Array, energy spectrum, mass composition, anisotropics 



1 Introduction 

Telescope Array (TA) is a cosmic ray observatory of the 
largest exposure in the northern hemisphere, located in 
the Millard County, 200 km southwest of Salt Lake City, 
Utah, USA. This is a hybrid experiment by collaborators 
from Japan, the US, Korea, Russia, and Belgium, using 
two types of detectors, surface detectors (SDs) and fluores- 
cence detectors (FDs). This can be regarded as a follow-up 
to the AG AS A and HiRes experiments, which used scin- 
tillation counters and fluorescence mirrors, to uncover the 
mysteries of ultra-high energy cosmic rays and to clarify 
their origin. In this paper, I review the major results from 
the TA first three years observation, including the energy 
spectrum, mass composition, and anisotropy of ultra-high 
energy cosmic rays. The activities for future extensions are 
also presented. 

2 TA detectors 

The TA-SD consists of 507 three square meters double- 
layered scintillation counters in a square grid with 1.2 
km separations in a total area of ~ 700km 2 (Figure []])• 
The average atmospheric depth at the TA ground level is 
~ 860g/cm 2 (~ 1400m above the sea level). Signals of 
shower particles are sampled at 50 MHz frequency and dig- 
itized with 12-bit FADC, not integrating all charges from 
photo-tubes in a ~ 10/iS time constant as in conventional 
air shower array experiments like AGASA. The SD data ac- 
quisition system is triggered by a coincidence, within 8/is 
window, of neibouring three SDs with signals greater than 



3 MIPs (minimum-ionizing particles) equivalent in their 
both layers. The SD array is fully efficient for cosmic rays 
with energies greater than 10 18,8 eV El. 

The fluorescence detectors are installed in the three sta- 
tions, those are the Black Rock (BR, 12 telescopes), Long 
Ridge (LR, 12 telescopes), and Middle Drum (MD, 14 tele- 
scopes) sites, viewing the sky above SDs. The 24 BR and 
LR telescopes were newly developed for TA, and the MD 
telescopes are refurbished HiRes- 1 detectors. One BR/LR 
telescope is comprised of a cluster of photo-tubes (or a 
"camera"), and a reflecting mirror of 3.3m diameter. The 
BR/LR camera has 16 x 16 2-inch hexagonal tubes cov- 
ering 15.5° x 18° of the sky. A telescope triggering sig- 
nal is generated when five adjacent tubes in a camera are 
fired by incident photons, and waveform data of all the 
cameras are recorded in a station [3 ]. The MD telescopes 
have 5.2m 2 mirrors and cameras of 256 pixels using 40mm 
Philips tubes. These telescopes are operated with the same 
electronics and data acquisition system of the HiRes- 1 ex- 
periment. The total FD fields of view are 108, 108 and 114 
degrees in azimuth and 3 ~ 31 degrees in elevation. See 
references d [H in this conference for details of the TA 
detectors. 

The TA detectors have been in operation since May 2008. 
The operational duty cycle of SD is ~ 98%, and the total 
exposure is ~ 2, 700 [km 2 str yr](~ 1.7AGASA) (~ April 
201 1). The telescope on-times are ~ 2400 hours (MD site), 
which gives ~ 1/3 HiRes- 1 exposure. The FD duty factor 
is ~ 11%, and it reduces to ~ 9% after good weather cut. 
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Figure 1: Layout of the TA detectors. 

3 Shower analysis 

3.1 Shower reconstruction with FD data 

The fluorescence detectors measure photons emitted from 
or around air showers with phototubes exposed to the night 
sky to determine the longitudinal profiles. The longitudinal 
profile of an air shower as a function of atmospheric depth 
X is written by the Gaisser-Hillas function, 
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the BR/LR data analysis, we have developed new codes 
suitable for the 10MHz sampling data acquisition system. 

Directly calculating the number of charged particles us- 
ing phototube signal and a given reconstructed air shower 
geometry is difficult due to the complexities of the FDs, 
as well as, accounting for the effects of Cherenkov light 
in the shower. For this reason, an Inverse Monte Carlo 
(IMC) method is employed. Using an IMC method, show- 
ers with varying X max and 7V max (or £ max ) are simulated 
and the number of measured photo-electrons is compared 
with that produced by the simulated showers. Simulation 
of the air shower allows for the treatment of Cherenkov 
radiation and ray-tracing methods are used to understand 
the acceptance of photons produced on the shower axis. 
Here we can include all the effects, such as photon shadow- 
ing by structures, mirror reflectances, and other calibration- 
related matters. We use independently developed FD simu- 
lation/reconstruction programs utilizing the inverse Monte 
Carlo method. From the differences between energies de- 
termined by the two codes, and together with ~ 5% and ~ 
3% uncertainties which originate from non-predetermined 
primary nuclear types and the missing energy correction, 
the energy uncertainty from the reconstruction procedures 
is evaluated as ~ 10%. Taking into account the system- 
atic uncertainties from FD calibrations [8] or atmospheric 
effects 0, the total uncertainty in the TA FD energies is 
- 22% Q. 

3.2 Measurement with the surface detector array 

Shower particles at the ground level of the TA site are sam- 
pled by the SDs with 3 m 2 area. From the lateral distri- 
bution of the shower particles determined from the local 
densities at SDsQ, an energy estimator, ^oo, the density at 
800m from the shower core is evaluated. The lateral distri- 
(1) bution function of shower particles we used is an AGASA- 
type formula ifTOlL 



where X max is the depth at the maximum development 
of an air shower, Xq is the first interaction point, A is a 
characteristic length of particle interaction. This equation 
was originally proposed to describe the development of the 
number of charged particles Q, and is also applicable to an 
energy deposit profile which can be directly measured with 
FDs. The calorimetric energy of the air shower is obtained 
by integrating the longitudinal profile, 
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where N max is the shower size at X max , (^) is the mean 
energy deposit of charged particles per unit length, and 
S(X) = £ m axf(X) is the energy deposit profile of the 
shower. The primary energy of the cosmic ray which in- 
duced the shower is obtained from E ca \ by taking into ac- 
count the missing energy which is an invisible component 
of the shower energy carried away by neutral particles. 

The data analysis procedures for the MD data are almost 
identical to those used in the HiRes experiment 0. For 



and the reconstruction procedures are adjusted to fit the 
TA SD data ifTTIl . In order to take into account differ- 
ent atmospheric attenuation of inclined showers, we con- 
struct a lookup table for conversion from the observables 
(£800, sec 0) to the energy E$r> by using a CORSIKA 
based Monte Carlo (MC) developed for TA ElElEl (Fig- 
ure [2]). The SD shower analysis program is tuned to repro- 
duce MC-thrown energies from the estimators (^oo , sec 0) 
of reconstructed showers. 

However, because of the lack of our knowledge of details 
of ultra-high energy hadronic interactions, there are rather 



1. The measured values are not necessarily proportional to 
number of particles: SDs measure mean energy deposit by shower 
particles, in terms of vertical muons. 



32nd International Cosmic Ray Conference, Beijing 201 1 




Figure 2: The SD energy chart 
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large systematic uncertainties in energy determination from 
the shower particle measurement at the ground compared 
to the calorimetric measurement by FDs. This is a long- 
standing problem in cosmic ray physics, and has been dis- 
cussed in comparisons between AGASA, an SD based ex- 
periment, and HiRes, which used FDs. In order to define a 
"unified" TA energy scale, we use hybrid events, which are 
detected by both SD and either of the three FD sites. From 
an Ef£ h - E^ h plot, where Ef* h and E^ h are ener- 
gies determined by the SD and FD reconstruction of the 
hybrid events respectively, we determined an overall scale 
factor < ^pD b /^D b >= 1/1-27 to obtain the energy of 
each SD event (Figure [3]). Therefore the energy of an SD 
event, which is not necessarily a hybrid event, is given by 
E =< E^ h /Ef^ h > E S B, where E S b is the energy de- 
termined from the SD reconstruction Q |2) . Note that the 
hybrid events of all the three combinations (SD and either 
of BR/LR/MD) are plotted in Figure [3] This indicates that 
the energy scales of all the three FD stations are in consis- 
tent even though the entirely different detectors and cali- 
bration schemes are used. 



4 Energy spectrum 

In Figure ffl I show three energy spectra above 10 18 eV 
from the TA data analyses, of the MD monocular mode 
ED, the BR/LR hybrid mode EH, and the SD analysis 0. 
Here the "BR/LR hybrid" means that we used FD data of 
BR or LR site, and also the shower particle arrival tim- 
ing at the position of an SD, which significantly improves 
the accuracy in determining shower geometries than in the 
monocular mode iflUl . For comparison, the AGASA lfl31l , 
HiRes [16] and the Pierre Auger Observatory (PAO) E3 
spectra are also displayed. We stress that all the three TA 
spectra are in good agreement (note that the SD energy has 
been scaled by 1/1.27 matching to the FD energy scale as 
described in the previous section). It is also important that 
the MD data, from the refurbished HiRes- 1 detectors, ana- 
lyzed with the identical procedures used in HiRes, is fully 
consistent with the HiRes- 1 and HiRes-2 results. 
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Figure 3: Scatter plot of the energies determined from SD 
and FD reconstructions using the TA hybrid events. Events 
of all the three combinations (SD-BR/LR/MD) are plotted. 

^[m-^tr-^-^V- 1 ] log 10 (£ a /eV) log 10 (£ c /eV) 



(2.4 ±0.1) x 10 



k 



18.69 ±0.3 



19.68 ±0.9 



-3.33 ±0.4 



-2.68 ±0.4 



-4.2 ±0.7 



Table 1 : TA-SD spectral parameters 

We applied a broken-power-law fit with two bending points 
to the SD spectrum, which is of the largest exposure in TA, 
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where k,£,m are spectral indices, E c are the bending 
points, and A is a normalization constant. The fit param- 
eters are listed in Table Q] The "ankle" structure is found 
at log 10 (^ a /eV) = 18.69, and the second bending point 
is also identified at \og 10 (E c /eV) = 19.68. The number 
of events observed above E c is iV" bs = 28, while expected 
^exp = 54.9 in an absence of the second bend. The signif- 
icance of the event deficit above the "cut-off energy E c is 
evaluated as Y%Lo Poisson(/i = 54.9; i) = 4.75 x 10" 1 ~ 
3.9a. 



5 Mass composition 

The information on the primary nuclear types of cosmic 
rays can be obtained from the longitudinal developments 
of air showers initiated by them. The fluorescence detec- 
tor is a powerful tool for this purpose, since it measures 
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Figure 4: TA energy spectra from SD (filled circle), FD 
mono (MD site, filled triangle), and BR/LR hybrid (filled 
square). 



shower profiles at different atmospheric depths with many 
photo-tube "pixels". Here we used BR/LR stereo events, 
which were detected at the two sites, because their shower 
geometries can be determined more accurately than in case 
of the monocular mode. We used X max , the maximum de- 
velopment point of an air shower, as the mass indicator: 
I max is larger (deeper in the atmosphere, lower altitudes) 
for protons or lighter nuclei, and smaller (higher altitudes) 
for heavier nuclei as iron. We compared the distribution 
of X max of observed showers and those predicted from a 
shower Monte-Carlo using the CORSIKA code G3. A 
caveat should be made here: since we observe cosmic rays 
at the fixed ground levels with the limited field of view, the 
distribution of observed X max should be different from a 
prediction by a shower generator like CORSIKA even if 
the primary composition is fully known, unless applying a 
sophisticated data selection @. We take another approach 
here. We use a "standard" data selection rule for the stereo 
events. Instead of finding a special selection filter for the 
data, we apply our detector and triggering simulation, and 
also the selection criteria to the CORSIKA showers exactly 
same as to the data fT9l . The distributions of X max of such 
Monte-Carlo showers with the detector simulation and the 
cuts can be compared with those of the real data. 

The distributions of X max of the observed showers for dif- 
ferent energy regions in E = io 18 - 4 ~ 19 - 6 eV are shown in 
Figure [5] [20]. The expected distributions for pure-proton 
and pure-iron primaries are also shown, obtained from our 
Monte-Carlo studies using the CORSIKA showers with the 
TA-FD detector simulations and the same data cuts . In or- 
der to examine the compatibilities of the observed X max 
and the predictions, KS tests were carried out calculat- 
ing the "p- values" as functions of energies. Here we used 
CORSIKA shower libraries of three different hadronic in- 
teraction models, QGSJET-01, QGSJET-II, and Sibyll. In 
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Figure 5: X max distributions of the observed showers of 
the different energy regions (markers), together with the ex- 
pected ones predicted from our Monte-Carlo studies using 
CORSIKA showers with the TA-FD detector simulations 
(histograms: solid for protons and dashed for iron nuclei). 
The QGSJET-II hadronic interaction model is used. 

all the energies and the hadronic interaction assumptions, 
the TA data are compatible with the proton primary predic- 
tions. In the higher energies (E > 10 19,4 eV), the p- values 
are also high for the iron hypothesis, but the statistics are 
still limited. 

I show the averaged X max as a function of energy in Fig- 
ure 01 20]. One can see again that the TA data is consistent 
with the proton-dominated composition of ultra-high en- 
ergy cosmic rays. Note that the proton and iron primary 
prediction lines in Figure are also obtained by applying 
the TA detector simulation and the same data cuts for the 
data, therefore both the data and the prediction lines in this 
figure cannot be directly compared with the results from 
other experiments, like HiRes or PAO. 



2. The PAO uses well-tuned data cuts which give a "true" X n 
distribution from the observed data |[T8l . 
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Figure 6: (X max ) plot as a function of energy. The proton 
and iron prediction lines are obtained through the TA-FD 
detector simulation and the cuts exactly same for the data. 

6 Anisotropics 



Figure 7: Auto-correlation analysis: the vertical axis is 
probability that the uniform distribution of cosmic rays ac- 
cidentally gives the same or a greater number of event pairs 
than the observed number of pairs in a given angular sepa- 
ration. 



The major difficulty in exploring the origin of cosmic rays 
is the magnetic deflections of charged particle trajectories 
from their sources to the Earth. However, anisotropies in 
the cosmic ray arrival direction distribution still give impor- 
tant clues, by finding hot spots of cosmic rays or correla- 
tions with possible source populations. We used the TA-SD 
events collected from May 2008 to April 2011, of zenith 
angles smaller than 45° and angular accuracies ~ 1.5°. 
Note that the SD event energies have been renormalized 
to the FD energy scale (1/1.27) as described in the Section 
13.21 The data set contains 854 events with energies greater 
than lOEeV, and 49, 20 for E > 40EeV, > 57EeV, respec- 
tively. 

First we examined auto-correlation of the TA events, or 
event clusterings as reported by the AGASA group ll2llL 
without making any assumptions on cosmic ray sources. 
The number of event pairs in a given angular window 5 is 
counted, and calculated a probability P(S) that the uniform 
distribution of cosmic rays yields the same or a greater 
number of pairs by chance than the observed one. The re- 
sult is shown in Figure [7] as a function of the window size 
ll22l . There is a rather small P(S) region around the an- 
gular scale c)~15 fori±?>57 EeV, but not statistically 
significant to support a deviation from the isotropy. 

Second, we searched for cross-correlation of the TA events 
and active galactic nuclei (AGN), which have been consid- 
ered as possible cosmic ray production sites. The proce- 
dures and parameters used in this analysis are similar to 
those used in the studies by PAO ll24l . We consider AGNs, 
QSOs and BL Lacs in the object catalog (the Veron-Cetty 
& Veron (VCV) catalog El) with redshift smaller than 
0.018, and search for correlation with the TA events of 
E > 57 EeV (20 events) within an angular separation of 
3.1° (Figure [8]). We assumed that all the objects have the 
same intrinsic cosmic ray intensities regardless of class or 
distance. The number of TA events correlating with the 



Event with the radius of 3.1 deg 




Figure 8: 20 TA events of E > 57 EeV and the AGN in the 
VCV catalog |23l . Galactic coordinates. 

VCV AGN is shown in Figure [9] as a function of the to- 
tal number of events with energies greater than 57 EeV, in 
an order of the detection date ll22l . We found 8 correlated 
events out of 20, while we expect 5 (25%) chance correla- 
tion if the source distribution were isotropic in the TA sky. 
The statistical significance of a deviation from the isotropy 
is also small in this analysis, therefore the AGN hypothesis 
on the cosmic ray origin is not supported by the present TA 
data. 

Lastly we investigated a correlation between the TA events 
and matter distribution in the universe. We employ the flux 
sampling method developed in |25l , and previously applied 
to the HiRes data (26). We use 109,408 objects from the 
2MASS Galaxy Redshift Catalog (XSCz) [27], which are 
located at distances 5 < d < 250 Mpc from the Earth, 
as tracers of the local matters in the large-scale structure 
(LSS) of the universe. The expected cosmic ray density 
map is constructed as a superposition of the contributions 
of the individual galaxies, taking into account an angular 
spreading (or smearing) as a result of the magnetic deflec- 
tions by the extra-galactic and Galactic irregular fields, flux 
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Figure 9: The growth of the number of TA events (E > 57 
EeV) correlating with AGN since 2008. 

attenuations due to cosmic ray interactions (e.g. photo-pion 
productions, e ± pair productions), and the TA exposure. 
Here the smearing effect is simply introduced as a Gaus- 
sian spreading with the smearing angle 6 S . The uniform 
component of cosmic rays, as cumulative contribution of 
sources beyond 250 Mpc, is also added. We assume that 
all the primary cosmic rays are protons, in accordance with 
the TA composition study 1 19, 20 ]. The obtained flux maps 
for E > 10, > 40, and > 57 EeV, together with the TA 
events, are shown in Figure Q0|. The compatibilities of the 
LSS hypothesis or the isotropy with the TA events were ex- 
amined with the KS test. The p-values as a function of the 
smearing angle are shown in FigureQj]for the structure hy- 
pothesis, as well as for the isotropic distribution l28ll29l . 
The TA data of E > 10 EeV is isotropic, and not compati- 
ble with the structure hypothesis. The p-values for E > 40 
and > 57 EeV are rather higher, but there are no significant 
differences between the structure and isotropy hypotheses. 
If we consider regular components of the Galactic mag- 
netic field, including a halo and a disk fields, the p-values of 
the structure model is mildly increased, while the compat- 
ibility with the isotropy is also maintained. Therefore we 
conclude that there are no convincing evidences of corre- 
lations with any classes of source candidates in the present 
TA data. 



Figure 10: The expected cosmic ray density map obtained 
from the algorithm in l25l with the smearing angle 9 S = 6° 
(contours). The dots are the TA events of the corresponding 
energies. Galactic coordinates (I = at the right edge of 
the figure, increasing t to the left, 180° at the center, and 
360° at the left edge). 
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7 Further developments 

There are several sources of energy uncertainties in the flu- 
orescence measurements, as atmospheric effects and fluo- 
rescence yields. There have been remarkable progress in 
laboratory measurements on fluorescence yields, however, 
it is not straightforward to apply the results to an actual cos- 
mic ray experiment taking into account all the atmospheric 
conditions. An in situ measurement of fluorescence light 
using known charged particle energies is of great advantage 
in a cosmic ray experiment, and the Electron Light Source 



Figure 1 1 : The KS test results for the compatibilities be- 
tween the hypotheses (structure and isotropy) and the TA 
data. No regular component of the Galactic magnetic field 
is considered. The horizontal line shows an a priori chosen 
confidence level of 95%. 
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Figure 14: TALE-FD field of view 



Figure 12: An example of the ELS beam shot seen from 
the TA fluorescence detectors. 
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Figure 13: Photo-tube signals obtained from an ELS trigger 
of the fluorescence detectors. 



(ELS) of TA is the unique one for this purpose Bill . This 
40 MeV electron shooter has been installed at the Black 
Rock site, 100m away in front of the fluorescence detectors, 
which enables us an end-to-end calibration of our detectors, 
atmospheric effects, and also fluorescence yields. The ELS 
has been operational since September 2010 (a beam shot 
example is shown in Figure IT2T). The signals (waveforms) 
obtained by the photo-tubes which see the electron beams 
from ELS are shown in \\3\ which are in good agreement 
with our Monte-Carlo expectations |3H . The analysis of 

+ U CT C A 9 • 2 4 6 8 J"* 12 A 14 + U 16 

the ELS data is now ongoing, m order to reduce the sys- 
tematic uncertainty in the FD energies by ~ 10%. 

I introduce $vo experimental actMtieJ^ndt 2 : def fclopments 
in TA. The one is the TA Low-energy Extension (TALE) 



. Precise ir^easuipmenjs of g energy s^ctnyji and mass 
composition of cosmic rays at energies 10 17 ~ 10 8 



are 



quite important by several reasons. In particular, a possi- 
ble transition of cosmic rays of Galactic sources to those 
of extra-galactic origins is of special interest. There is 
a general agreement that heavier nuclei are dominant in 
cosmic rays around the knee region {E = 10 15 ~ 16 eV), 
and the average mass increases with energy up to ~ 10 17 
eV. This is consistent with the supernova-origin scenario 
of cosmic ray production and acceleration in our Galaxy. 
On the other hand, there are several experimental results 
which suggest proton dominant composition above 10 18 
eV |[33l . Therefore it is expected that there could be a dras- 
tic change in mass composition of cosmic rays at 10 17 ~ 18 
eV, which gives a crucial experimental constraint on astro- 
physical interpretations of the proposed second knee or the 
ankle structures in the energy spectrum. TALE is designed 
to lower the minimum detection energy of TA by about an 
order of magnitude, with the hybrid detection technique. 
Reconditioned HiRes-2 detectors as the TALE FD will be 
installed in the immediate vicinity of the MD site i.e. the 
HiRes-1 detectors. The TA-TALE FD covers from 3° to 
59° in elevation (Figure [14]), which enable us to detect 
lower energy cosmic rays and a bias-free X max measure- 
ment. An infill array is being installed in front of the TA- 
TALE FDs. A sub-array of 45 counters with a spacing of 
400m at a distance of 1.5 km from the TALE-FD gives a 
10% efficiency for cosmic rays at 3 x 10 16 eV. In order to 
cover more area, we change the spacing to 600 in a range 
of 3 ~ 5 km from the FD. For larger distances we use the 
spacing of the main TA SD array, 1200 m. The layout of 
the TALE SD is shown in Figure [T31 

The another experimental extension in TA is a bistatic 
radar detection of ultra-high energy cosmic rays (TARA) 
Ol [351, firstly explored by the MARIACHI project 1551 . 
This is a lower-cost, 24 hour duty cycle remote sensing 
technique for the next generation cosmic ray observato- 
ries. An important advantage of the bistatic radar tech- 
nique among the several radio detection methods proposed 
for cosmic rays comes from the fact that the scattering 
cross-section is largest in the forward direction [37]: a ra- 
dio pulse sent from the transmitter station, and reflected 
by air shower plasma can be effectively detected at the re- 
ceiver station. The signal of the reflected radio pulse has a 
"chirp" waveform due to a Doppler-like effect by the scat- 
tered which is the unique feature of the air shower echoes. 
We have installed a 2 kW analog television transmitter do- 
nated from KUTV Channel 2 in Salt Lake City (FigureQS]). 
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Figure 15: Layout of the TALE infill array. The light- 
colored circles indicated with "A"-H" or 4 digits + "t" are 
the TALE counters to be added |[32l . The other darker cir- 
cles with 4 digits + "T" represent the existing counters in 
the north-west part of the TA SD H). 

This was commissioned in January 201 1, operated at a sin- 
gle carrier frequency at 54.1 MHz. The receiver antennas 
have been deployed at the TA-LR site, 50 away from the 
transmitter station (Figure Q21Q2]). We also plan to utilize 
the radar transmitter in conjunction with the TA ELS. The 
transmitter will be upgraded at least to 40 kW by the sup- 
port from the US-NSF/MRI. 

8 Discussion 

There is a ~ 20% energy scale difference between TA and 
PAO, as inferred from the energy spectra. It should be 
noted that there are 22 ~ 23% systematic uncertainties in 
energy determinations in both experiments, therefore it is 




Figure 16: TA bistatic radar (TARA): the radio transmitter 
and the receiver have been installed at a distance of 50 km. 



Figure 17: TA radar receiver 

possible to say that the two results are consistent within 
the systematic errors. However, the 20% difference is not 
negligible and cannot be overlooked for better understand- 
ing of air shower experiment and astrophysics of cosmic 
rays. It is partly attributed to the different fluorescence 
yield models employed in the two groups (~ 10%), but not 
fully explained only with it. Some cross-calibration of the 
detectors will be required to reduce the energy scale differ- 
ence, e.g. exchanging the calibration tools or using the TA 
ELS (this is the author's personal opinion, and it has not 
been officially discussed yet between the group headquar- 
ters). 

Both the TA and PAO (X max ) results show that proton 
is dominant at energies above 10 18 2 eV. It can be inter- 
preted that we have not observed an increase of average 
masse of cosmic rays in the lower energies, which is ex- 
pected from the results of the air shower array experiments 
targeted at and above the knee region. This indicates the 
importance of the lower energy extensions like TA-TALE 
El, or AMIGA and HEAT in PAO (38|. It is controver- 
sial that there is an apparent difference in the composition 
studies in TA and PAO in higher energies. The TA (X max ) 
result is consistent with a proton-dominant composition of 
cosmic rays |[20l (there is one event with a rather smaller 
(I max ) at the highest energy bin, but it is of essentially no 
significance), but an increase of average cosmic ray mass 
has been suggested from the PAO data (e.g. CU[39l|40|). 1 
propose that we try the same analysis procedures: now we 
are developing a bias-free data cut and determining (X max ) 
as carried out in PAO (our preliminary result shows that this 
brings same amount of shifts in (X max ) to both the data and 
Monte-Carlo, therefore the conclusion of proton-dominant 
composition is not changed at all.). It is possible to em- 
ploy the TA-style analysis given in (191 [20| in PAO. It is 
also important to revisit the high-energy hadronic interac- 
tion models, after the recent experimental progress by LHC 

ED. 
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In 2007, the PAO claimed that the distribution of cosmic 
ray arrival directions is not isotropic and suggested a pos- 
sible correlation with AGN B2l . However, the statistical 
significance of the AGN hypothesis in the PAO data has 
been decreased in their updated analysis O |40). The 
current TA data shows no statistically significant small or 
large-scale anisotropies. Unfortunately the interesting re- 
gion around Cen A is out of the TA sky. There is no hot 
spot around the Virgo cluster in the TA data. The absence 
of apparent anisotropy or correlation with known astro- 
nomical objects could bring a difficulty in the interpreta- 
tion of the spectral features, in particular of the steepen- 
ing at ~ 5 x 10 19 eV, firstly claimed by the HiRes ex- 
periment El, and also found in PAO [45] and TA 0. If 
the primary cosmic rays are protons, and the steepening is 
due to cosmic ray energy losses through interaction with 
the cosmic microwave background photons (the Greisen- 
Zatsepin-Kuzmin effect l46l ). If this is the case there could 
exist a cosmic ray "horizon" at a certain distance from the 
Earth, and hence a strong correlation with nearby objects. 
Another interpretation is that a heavier component is domi- 
nant at the highest energies, therefore the steepening we ob- 
served is not due to the GZK effect, but to other processes, 
e.g. acceleration limit at production sites. In any case, no 
consistent picture of cosmic ray nature can be drawn from 
the current observation data, although we believe that we 
are just about the answer. Further efforts in data analysis 
are required in determining energies and composition, and 
statistics for anisotropy studies. 

9 Conclusions 

Telescope Array is the largest cosmic ray observatory in 
the northern hemisphere, which has been operational in 
full configuration since May 2008. The TA exposure is 
~ 1.7 AGASA in SD, and - 1/3 HiRes in FD. We have 
developed detector simulators and shower reconstruction 
programs for both SD and FD, and the distributions of the 
observed data are in excellent agreement with expectations 
using a previously measured spectrum of cosmic rays. The 
energy scales of the three FD stations (BR, LR and MD 
sites) are also in agreement, even though different detec- 
tors and calibration schemes are used. 

I presented the major results from the recent progress of 
the TA science, including the energy spectrum, mass com- 
position, and anisotropy of ultra-high energy cosmic rays. 
The energy spectrum of cosmic rays has been derived from 
analyses of the FD, SD, and the hybrid data. All of the three 
spectra are in good agreement. From the spectral fit of the 
SD spectrum, which is of the largest exposure in TA, we 
found the ankle structure at 10 18-69 eV. We also identified a 
steepening in the spectrum at 10 19,68 eV with a change of 
the spectral index from —2.68 to —4.2. An extended spec- 
trum of cosmic rays at the highest energies is ruled out by 
a statistical significance of 3.9a. We studied the distribu- 
tion of X max , at which an air shower has its maximum size, 
using the FD data. We compared to the X max distribution 



of the data and expectations from proton or iron nuclei pri- 
maries. The TA data is consistent with the proton-dominant 
composition. Furthermore, we examined the arrival direc- 
tions of cosmic rays with energies greater than 10 EeV, and 
searched for anisotropies or correlations with known astro- 
nomical objects. We concluded that no statistically signifi- 
cant anisotropies are found in the present TA data. I intro- 
duced a challenge for an end-to-end calibration of the FDs 
using the electron accelerator (ELS). I also discussed the 
further experimental developments in near future, the TA 
low energy extension (TALE), and the bistatic radar detec- 
tion of ultra-high energy cosmic rays with TARA. 

It will have been 100 years in 2012 since the discovery of 
cosmic rays by Victor Hess. The world financial and po- 
litical clock is ticking. A critical question to cosmic ray 
physics is that : "who knows we would discover the origin 
of cosmic rays?" If the TA energy measurement and the 
composition study are correct, the final unsolved problem 
is anisotropy. Our study shows that the statistical power of 
distinguishing different hypothesis on cosmic ray origins 
are still weak, and smaller than 50% in the current TA data 
|[28l . This will be significantly improved if we have statis- 
tics doubled or tripled, therefore we need observations sev- 
eral more years. We'll also have to discuss the next gen- 
eration cosmic ray observatory with a huge coverage, to 
explore the astronomy by charged particles. 
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